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A CARDIOVASCULAR SYSTEM MODEL FOR 
LOWER-BODY NEGATIVE PRESSURE RESPONSE 

Baker A. Mitchell, Jr. 

Robert P . Giese 
DEPARTMENT OF BIOMATHEMATICS 
THE UNIVERSITY OF TEXAS 

M. D. ANDERSON HOSPITAL AND TUMOR INSTITUTE AT HOUSTON 
1. INTRODUCTION 

One method for studying complex physiological control systems 
is through the use of mathematical modeling. This approach permits 
investigators to identify gaps in present knowledge of a subject as 
well as permitting a basis upon which to place new findings . The 
ultimate goal of these simulation systems is to provide a source of 
predictions for those situations when costly and timely research is 
prohibitive. Therefore, this type of effort finds wide application 
in dealing with physiological problems and space flight. 

The purpose of this effort was to modify, for lower-body negative 
pressure studies, a mathematical model of the cardiovascular system 
which this contractor has previously developed. Since the mathe- 
matical model was about twice as large as the capacity of the EAI 680 
computer available, an all-digital computer model was written for 
installation on the Univac 1108. This all-digital computer model was 
specifically written to allow orderly, straightforward expansion to 
include exercise, metabolism (thermal stress), respiration, and other 
body functions, the accuracy of which are highly dependent on a good 
cardiovascular fluid dynamics model. 
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2 . TERMINOLOGY 


Since the word Model is so ambiguous, an explanation of its 
usage in this paper is in order. 



The Mathematical 
Formulation of 
this idea 



o- 


H=hzh=H 


The Computer 
Approximation of 
this Formulation 



The human cardiovascular 
system under all possible 
environments . 


The human cardiovascular 
system under LBNP , gravity, 
tilt angles, deconditioning, 
etc. 


The Abstract Model 


V - V /(F in - F .»t> dt 

The Mathematical Model 


DO 481 I - 1,15... 

The Computer Model 


The Computer 
Approximation 
with a given set| 
of parameters 



Vol Head = 500mi. 

A Specific Computer Model 


A Particular 
Experiment performed 
with these parameters 



An Experiment or Computer Run 


Figure la 


Figure lb 
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Figure la represents the modeling process in general. Each 
step is a further abstraction of the step above. Figure lb represents 
the corresponding steps in the cardiovascular modeling project. When- 
ever the word model is used in a place that could be ambiguous, the 
appropriate modifier abstract , mathematical , computer or specific 
computer will be used. 

3. THE HUMAN CARDIOVASCULAR SYSTEM 

We cannot hope to model every function of the human cardio- 
vascular system. We can, however, choose a sufficiently large 
subset of its functions so that certain experiments on the model 
will be predictive of the behavior observed when the same experiment 
is performed on the human being. 

4 . the part of the cardiovascular system modeled 

Only the relatively short term changes of the cardiovascular 
system are to be modeled. There are longer time constants not 
accounted for, such as: 

A. Vena Cava stress relaxation with time constant in the 
order of 12 minutes. 

B. Intravascular and extravascular fluid shifts with time 
constants of the order of 90 minutes. 
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C. Metabolic and biochemical stress with various time constants. 

5. THE ABSTRACT MODEL 

The abstract model configuration is shown in figure 2. It attempts 
to match the behavior of the fluid dynamics of the following anatomical 
elements. The number associated with each element is the same as is 
used in the computer model to follow. 

A. Fluid Compartments 

1. Heart Chambers 

8. Right Ventricle 

11. Left Atrium 

12. Left Ventricle 

2 • Aorta 

13. Ascending 

14 . Carotid 

15. Thoracic 

3* Other Fluid Compartments 

1. Head 

2 . Arms 

3 . Trunk 

4. Kidney 

5 . Gut 

6 . Legs 

7 . Vena Cava 
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9 . Right Lung 
10. Left Lung 

B . Other Active Elements 

1. Carotid Baroreceptor 

2. Cardiac Pacemaker 
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Figure 2 
ABSTRACT MODEL 



NOTE : 


Large numbers are 
Small numbers are 


organs ; 
flow paths 
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This abstract model considers the cardiovascular system to 
consist of 15 separate compartments with the 21 associated flow paths 
interconnecting them as in figure 2. All of the body or model fluid 
is considered to be in these 15 compartments, and no blood loss 
(bleeding) or intra-vascular or extra-vascular fluid shifts are 
considered. Since respiration and energy usage. are not explicit 
in this model, there is no oxygen and CO^ content of the blood 
considered. Thus, heart rate for the model is entirely determined 
by pressure in the second aortic chamber (baroreceptor) , and not 
secondarily by O 2 and CO 2 concentration at chemoreceptor sites as 
is the case in the real world. 

Since only one pressure is associated with each compartment, 
some simplifications are necessary. For example, while one pressure 
is sufficient in the description of a heart chamber or a section of 
the aorta, it is not totally sufficient in the description of the 
legs, head, etc. Thus, the pressure computed is that in the 
arterioles supplying those organs. This convention will be especially 
useful for future secondary control of heart rate by 0 2 and C0 2 and 
for regulating of blood flow by the various organs. 


6. THE MATHEMATICAL MODEL 


The equations which describe certain aspects of the behavior 
of each compartment of the abstract model will be described along 
with the behavior characteristics pertinent to the simulation. 
Symbols, units, and expressions for physiological constants are 
all contained in the appendix. 
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A. 


General Fluid Compartments 


flow 

Volume , 

flow 

P 

flow 

p 

in 

Pressure 

out 

m 


n 


Figure 3a 


Figure 3b 


Figure 3 illustrates the basic equations for flow, 
pressure and volume. 


1 ) 


v 



F ) dt 
out 


* 


2 ) 

3) 



If gravity is considered these become: 


4) 

5) 



R 

m , n 


*See Appendix v for a description of these symbols. 
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Equation 5) can be further refined to consider the 
inertia of the blood by using the following relationship 
from Mitchell, converted here to our notation: 


6 ) 



R * 
m, n 


F + L * 


dF 

dt 


which becomes 
after adding 
gravity , 


dF 

P - P - G = R * F + L * 
m n m , n dt 


dF 

F * R =P - P - G- L — 
m , n m n dt 


or. 


7) F - 


dF 

P m - P n - 6 - L ^ 


m , n 


Thus, equations 1), 4), and 7) become the basic 
equations of the model. Equation 1) maintains fluid 
mass conservation in the compartments. Equation 4) 
describes the pressure - volume relationship of the 
tissue mass forming the organ, and equation 7) is the 
classical pressure-flow relationship considering gravity 
and inertia. 


B. Heart Chambers 

Heart chambers are described by equations 1) , 4) , 
and 7), with C and R as variables. Varying C as a function 
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of time to simulate "pumping" has become common practice 
and seems to yield qualitatively satisfactory results. 

8) C = C(t) 

Resistance of a valve is determined by the direction 
of flow through the valve . 


9) 


R(F) 



00 


F >> 0 

for perfect valve. 

F < 0 


All heart chambers are governed by equations 1) , 4) , 
7) , 8) , and 9) , with the appropriate values of R and C(t) 
for each chamber being determined empirically. 


Aorta 

Beginning with the stress-strain relationship for a 
thin-wall cylinder: 

10) P- = eAi , 

r r 

the following equation governing incompressible fluid flow 
in elastic tubes may be derived: 

Ed dA 

11) P = E ( A - A ) + G + — , 

sO A dt 

where E is static elasticity and Ed is dynamic elasticity, 
s 

The derivation of their equation is found in appendix I . 
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Note that 11) is simply 4) with the addition of 
dynamic term (Ed/A) * dA/dt. 

Equation 1) for conservation of mass and equation 7) 
for pressure versus flow together with 11) complete the 
description of the aorta. 

The aorta was discretized into three 10cm segments 
and the constants in equations 1) and 7) were calculated 
accordingly. 

Carotid Baroreceptor 

The vagus nerve impulse rate is related to aortic wall 
strain which may be obtained from equation 11) . The exact 
nature of the transfer function between wall strain and 
impulse rate is presently under investigation. 

Cardiac Pacemaker 

Again, exact relationships which are physiologically 
justifiable are presently lacking. Here again "empirical 
curve fitting" is a more apt description than is the word 
simulation . 

External Forces 

Hydrostatic pressures are added to all fluid compart- 
ments; the heights of the pressure heads are based on the 
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body being in a straight posture. Additional acceleration 
terms may be easily added. 


7. THE COMPUTER MODEL 

A. Program Design 

The LBNP Program is structured in levels of activities, 

most of which are nested as shown in figure 3. A computer 

job is a complete unit and may contain one or more runs. Each 

run is self contained, but is broken down into one or more 

time segments. The time segments represent intervals of 

time during which the model is using a fixed set of parameters. 

Body parameters such as resistance, LBNP, compliance, etc., 

I may only be changed between time segments. Thus, to simulate 

an experiment, several time segments, or even several runs, 

may be required. Each time segment contains one or more 

buffers of 200 sets of output variables. For example, if 

100 sets of output variables per second were desired, then 

the length of a buffer would correspond to two seconds of 

model time. Each time segment must then be an even number 

t 

of seconds, and there will be - buffers in a time segment 
of length t. As mentioned above, there are always exactly 
200 output units in a buffer. Each output unit contains one 
or more compute units. In general, the stepsize At between 
compute units is considerably smaller than the resolution of 
^ the output device. 
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Compute Unit 
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98 Card 






The above levels of activity are all nested as shown in 
figure 3. There is one activity that does not nest with these 
other units; i.e., the cardiac cycle. The heart begins to 
beat at time=0 with the beginning of atrial contraction. 

Since there is essentially a continuum of heart rates, and 
heart rate may change after each beat, there is no assurance 
that the beginning of any heart beat except the first will 
coincide with the beginning of any other time segment, buffer, 
output unit, or even a compute unit. Thus, the cardiac 
cycle can not nest with the other units since it is not 
properly contained in any other unit, nor does it properly 
contain any other unit. * 

B. Subcomponent Detail 

For clarity, the subcomponent description will begin with 
the innermost level of activity and work outward. 

1. The Compute Cycle 

The compute cycle is controlled by an internal 
"clock" which runs continuously for each submitted 
job. This clock is represented internally by several 
variables within the program to insure: 1) no float- 

ing point found off or "drift"; 2) no integer over- 
flow or "involuntary resetting"; 3) accurate sub- 
division as small time units to insure accurate 
integration; 4) separate timing of each heart beat 

*See Appendix III for flowcharts. 
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so that no parameter change can occur in the middle 
of a heart beat. This "clock" may seem more complex 
than necessary, but each part is necessary to insure 
1) and 4) . Both integer and real variables are used, 
and the basic units of the clock are seconds and 
milliseconds. In addition, many variables that are 
used often in the program such as ^ , ^*A t ’ At ’ 2 * At> 
are pre-computed to minimize computation time and to 
use multiplication whenever possible instead of 
division. In-line comments in the program define 
both the working of the "clock" and the units of 
these associated variables. 

The compute cycle as flowcharted in Appendix III 
uses the equations : 


1) V = V Q + /(F ln - F out )dt 

2) p = c ~ P r + G » 


and 


7) 


dF 



m , n 


along with the associated equation of the heart and 
aorta stated previously, to compute new valves for 
each V, P, and F, from the values of the last compute 
step. * 

*Derivation of Aortic pressure volume equations is 
found in Appendix I. 
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In equation 4) the value of C is a function of 
time for the heart chambers to simulate pumping. For 
the other organs, C is actually a function of volume; 
although it is often considered a constant for small 
volume changes. Computation speed was increased 
considerably by writing a FORTRAN function that 
returns a value of P for a given volume. In the 
function C, the value of P is computed by linear 
interpolation of a table of five pairs of values of 
F and V. 

Equation 7) is restated as follows: 



n n-1 

Substituting: F F for d_F 

At dt 


n n-1 
F - F 



Rearranging: 


n 

R F = P 
m , n 


m 


- P 


n 


G 


L 

At 


n 

F + 



5 
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R 

m , n 


L n 

+ aT f 


p 

m 


- P 


n 


n-1 


G + 


L F 
A t 


» 


n 

F 


P - P - G + 
m n 


n-1 
L F 

At 


( R + — ) 
m,n T At 


This last formula turned out to be acceptable 
in actual computation. 

All variables are computed in standard medical 
units. Several subroutines are included with the 
program to convert between scientific and medical 
units if this becomes desirable. A symbol table and 
list of included functions and subroutines are 
included with this report. In-line comments in both 
the calling program and the subprogram explain the 
use of each subprogram. 

Valves were placed in the same flow paths as in 
the human body. A valve simply prevents reverse flow 
through a flow path. The valves in this program are 
"perfect" in the sense that they allow no back flow 
in closing. Similiarly, all volumes are constrained 
to be positive. If the computed flow is found to 
drive a volume negative, the volume is set to zero 
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and the flow recomputed. 


Stepsize can be a critical factor in this 
program, especially if relatively unrealistic body 
constants are used. In general, 200 steps per 
second are sufficient, and in most cases the results 
are indistinguishable from a higher compute frequency. 
In varying only stepsize, the most striking difference 
between 200 and 3200 steps per second is that in the 
latter case the compliance function of the heart is 
more nearly correct. This yields a slightly more 
efficient heart, resulting in about four beats per 
minute decrease in heart rate. If a stepsize is too 
large because of a compute frequency of 100 or less, 
then a cyclic ringing often occurs in the aorta. 

This will be especially evident with larger than 
realistic flow rates in the aorta, very high cardio- 
vascular output, or larger than real-life inertia 
term for the blood flow. To some extent this ringing 
can be artificially dampened by introducing a more 
elastic aorta using a different compliance curve 
or by increasing the inertia term on the aorta wall. 
However, using a higher compute frequency, and thus 
smaller stepsize, is probably the most accurate way 
of handling this difficulty. Note, however, that 
over 90% of the executive time of this program is 
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within this loop, so doubling compute frequency 
nearly doubles running time. Thus, the indiscriminate 
use of a very high compute frequency to eliminate any 
chance of this ringing will use up large amounts of 
computer time. 

2. The Output Cycle 

Since the desired compute frequency may not be 
the same as the desired display frequency, the program 
allows the user to specify an output frequency 
different from the chosen compute frequency. 

To simplify the program we require that the 
output frequency be an integer power to two times the 
buffer size of 200. Thus, acceptable output 

frequencies are , 50, 100, 200, 400, , points 

per second. Output frequency should always be 
adjusted considering the resolution of the output 
device. Compute frequency must always be a non- 
negative, power-of-two multiple of the output frequency. 
Thus, for a frequency of 100, a compute frequency 

should be 100, 200, 400, . In many cases it is 

desirable to have a relatively small output frequency 
(such as 50) and a relatively large compute frequency 
(such as 400) . 
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3. The Output Buffer 


All 51 output values of pressure, volume, and 
flow terms are stored in a 51 x 200 output buffer, 
and are output to magnetic tape as a single unit. 
This allows for: 1) separation of compute logic 

from graphical output logic so that graphical output 
units may be changed at any time with a minimum of 
effort; 2) the display of as many or as few output 
parameters as is desired; 3) re-play of a job 
utilizing different scaling factors or outputting 
different variables without the overhead of re- 
running the entire compute program; 4) minimize 
core requirement of the program without using an 
overlay structure. 

The current value of each of the 10 derived 
variables are also written out at the end of this 
51 x 200 buffer. The computation of these variables 
is discussed in the following section. 

After each buffer is written to tape, a check 
for constant blood volume is performed. This check 
is made for accumulated truncation and round-off 
errors; and the excess or shortage of blood in the 
entire system is distributed over all 15 organs in 
the same proportion as the original blood volume in 
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those organs. For example, if there is a lm£ 
excess in a 5000m£ system where the arms originally 
contained 500m£ of blood, then the arms would loose 
0.1m£of excess blood. The total blood volume before 
each update is printed on the computer listing. In 
most cases, the accumulated error at time of update 
is about 1 or 2m£, or considerably less than 0.1%. 

4. Cardiac Cycle 

Circulation is produced in the model in the 

following manner. For the pressure formula of the 

V 

heart chamber, P = — - P + G, the compliance, C, 

Li i 

is a function of time. The following flow diagram 
illustrates the logic that determines the heart rate 
(figure 4 ) . 

Note that in all cases the heart is allowed to 
complete its current beat before starting the next 
beat. The maximum rate allowed with this model is 
180 beats per minute, since above that point the heart 
does not fill rapidly enough and stroke volume is 
decreased to the point that cardiovascular output 
drops beyond that point. 

No minimum heart rate is imposed on the model, 
except that no negative rate is allowed. Under some 
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Figure 4 







gravity conditions with certain body angles the 
heart can be made to stop. This is physiologically 
as correct as is possible without carrying oxygen 
and CO 2 content to the blood. 

5. Derived Variables 

Systolic, diastolic, and mean aortic pressure 
are computed in the following manner. The heart rate 
is converted to number of compute steps per beat. 

Then the heart is allowed to complete a beat with 
that number of compute steps. For example, if we have 
200 steps per second compute rate and a heart rate 
of 60 beats per minute, this one beat is 200 steps. 
Pressures are computed as follows: 

n 

Systolic pressure = max (Pressure at step i) ; 

i= 1 

n 

Average pressure = Z (Pressure at step i) ; 

i = l 

n 

n 

Diastolic pressure = min (Pressure at step i); 

i=l 


where pressures are measured in the second 10cm of 
the aorta ( P-^ in the notation of the computer 

model) . * 

*See Appendix II for description of computer variables. 
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6. Time Segments 


A time segment is the basic interval of the 
program used for specifying a sequence of events in 
simulating an experiment. Since there is no require- 
ment for this all-digital program to be run real- 
time as in the case of most analog or hybrid models, 
a system was designed to allow the program user to 
design an experiment in advance of the computer run, 
in much the same way as a programmer lays out the 
logic of his computer program. 

Each time segment is preceded by the reading of 
any number of parameter cards to define a physiological 
personality and the external forces to be applied to 
the model. When a "time card" 

( 95 20 . 

such as 1 the above example is read, the time segment 
will start and will run for the designated number of 
seconds. At the end of that time segment, new 
parameter cards are read until another "time card" 
is read. 

Since a parameter change can affect the heart 
compliance, there is a delay factor in the time 
segment loop that will not use the new parameters 
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until the beginning of the first complete heart 
beat within the time segment. Note that each time 
segment is exactly the length specified on the "time 
card", and no error accumulates because of this delay. 

7 . Runs 

Each run may simulate an individual experiment 
or portion of an experiment. The results of the 
simulation are then stored on an intermediate magnetic 
tape, and can then be displayed at any time without 
requiring the entire simulation to be rerun each time. 
An example of the use of this feature follows. The 
user wishes to vary the compliance of one organ 
(without loss of generality, say the legs) and observe 
the effects on heart rate and flow to the brain.* He 
may then display only a few desired variables such 
as flow to brain, heart rate, and leg volume. He 
then repeats the simulated experiment several times 

using the same deck setup until the desired behavior 

\ 

is observed. The user may then want to display all 
variables requiring several passes from the same 
intermediate tape without rerunning the entire 
simulation. Finally, if he is satisfied with the 
results, he may again play back the intermediate tape 
and display the parameters at a scale suitable for 
publication. In this way, the investigator may 

*See Appendix IV for a description of compliance cards. 
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display at \ sec or 1 sec per inch to show fine 
details of a wave form, or compress long periods 
of time into a short space by using a display of 
20, 60, or even 600 seconds per inch. 

The deck setup for both the experiment and the 
display is explained in the User's Guide which is 
enclosed with this report. Full multi-job capability 
is included in the program for the user who wishes 
to run several experiments with the same computer 
job . 

C. Input-Output 

1. LBNP Program 

a. Input 

Control cards and deck setup are 
discussed in the User's Guide.* The entire 
set of cards that generated the output sample 
included in the user's guide is also attached 
to that document. All initial pressures, 
volumes, compliance, heart rate, resistance, 
inertia values, body angle, gravity, compute 
and output frequency, and body dimensions 
must be input to the program. In general, 
there are no default options for these 
parameters. However, since data cards may 

*See Appendix VI for User's Guide. 
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be read into the LBNP Program in any order, 
and since duplicate cards may be read 
retaining only the value of the last card, 
there is an easy way to avoid such omissions. 
One needs only to input the entire data deck 
supplied with the program and add modifica- 
tions at the end of that deck. For example, 
if more conditioned legs were desired, a 
second card 6 could be included behind the 
enclosed deck. 

b. Output 

All input data constants are output at 
the beginning of each run. Run number, time 
segment number, and length of time segment 
in seconds is output for each time segment. 
Heart rate, cardiovascular output, systolic, 
mean, and diastolic pressures are output for 
each heart beat. For the first time segment 
of each run, the flow, volume, and pressures 
11 through 15 are output in their entirety. 
This is to allow the user to see accurately 
computed values of this critical part of the 
cardiovascular system without waiting for the 
plot to be returned. 
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After each output buffer is written 
out on tape, the following information is 
output: first line, segment number and 

total accumulated time within run; next 21 
lines, flow, pressure, and volume of the 
model at the last compute step; last line, 
total volume in system before automatic 
volume update. 

2 . DISPLAY Program 
a. Input 

The first 62 cards of this deck must 
be included in the order given for each run. 
These cards labeled starting with 101 to 411 
are scaling factors for all output variables. 
The hundred's digit of 1 is for flow, 2 is 
for volume, 3 is for pressure, and 4 is for 
special output parameters as discussed in 
the user's guide. The field 5-8 is for 
starting value of the graph, and field 9-12 
is scaling factor of data to plot dimension. 
Examples of usage are given on the following 
page . 
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Field 

5 6 7 8 9 X 0 

Minimum 

Value 

Maximum 

Value 

0 

. 0 

1 

0 

100 

0 

. 0 0 

5 

0 

200 

0 

- . 0 

1 

0 

-100 

5 0 

. 0 

1 

50 

150 

9 0 

. 0 

5 

90 

110 

Minimum 

Value 

is 

field (5-8) 


Maximum 

Value 

is 

field (5-8) 

1 

+ field(9-12) 


For each job on the intermediate tape 
to be read by this program, the user must 
submit cards to indicate what fields are to 
be displayed, and which vertical and horizontal 
scaling factors to use. The user's guide 
explains the necessary format and usage. This 
display program may be used on either a XDS 
Sigma 5 computer or a Univac 1108 computer, 
with either a 12 or a 22 inch plotter. 
Alteration between these computers and plotters 
is accomplished by adding a c in column 1 of 
some program cards , and removing a C in column 
one of others. In-line comments locate those 
areas and explains the needed changes. There 
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should be little or no difficulty changing 
this program to other computers or other 
size plotters. 


b . Output 

The complete printer output may be 
ignored since it contains information of 
little use to the program user. The output 
contains only a listing of the values on the 
scaling cards, plus the value of the para- 
meters G(l) through G(10) for each output 
buffer. These variables are: 


G ( 1 ) IS L.BNP 

GO) IS SYST8LIC PRESSURE 
G ( 3 ) IS DIASTBLIC PRESSURE 
5(4) IS' MEAN BL09D PRESSURE 
G( 5 ) IS HEART RATE 
G ( 6 ) IS B0DY ANGLE 
G ( 7 ) IS CARDIBVASCULAR 9UTPUT 
G ( 8 ) IS GRAVITY 

3(9) T0 G ( 10 > ARE CURRENTLY N9T USED. 


An example of the plotter output is 
included in the enclosed user's guide. Up 
to 20 graphs may be displayed per computer 
run (up to 40 with a 22 inch plotter) . The 
graphs may be displayed in any order and the 
time scale is continuously variable. The 
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graphical display is especially easy to 
read if it is plotted on graph paper 
divided either 10 divisions per inch or 
10 divisions per half inch. 

THE SPECIFIC COMPUTER MODEL 

The set of input data constants and functions that compose the 
specific computer model are realistic for an 80kg man 180cm in height, 
in good physicial condition. Total blood volume is just under 5h 
liters and cardiovascular output prone and at rest is about 80m£/sec, 
or about 5 liters per minute. A resting heart rate is about 72, and 
systolic, mean, and diastolic pressures are approximately 140, 100, 
and 80 respectively. 

The distribution of blood in the model is approximately correct; 
however, it should be calculated more accurately if very close 
comparison to a particular human subject is desirable. This fixed 
model responds to G loading and lower body negative pressure in much 
the same way as a human would for short periods of time . Longer range 
changes may be simulated by a change of compliance, especially in the 
legs and vena cava. The exact parameters for this model are given in 
the user's guide. 

AN EXPERIMENT 

Although a very large number of experiments have been run with 
this program, only a single experiment is included with this paper. 
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The experiment involves applying lower body negative pressure to the 
prone fixed computer model of section 8. Sixteen seconds were allowed 
for the model to stabilize with zero negative lower body pressure (LBNP) . 


Then - 30 , - 40, - 50 mm of mercury were applied for 8 seconds each, 
followed by a 16 second run at zero LBNP. The deck setup as well as 
the graphical output is included with the User's Guide. 

10. EXAMPLES OF PARAMETER CHANGES * 

The user's guide included with this report contains a listing of 
parameter cards used for the fixed computer model. In this section we 
illustrate by example various changes in the model. 

A. More Conditioned Legs 

With better conditioning, the legs are less compliant 
and will accept less blood upon standing or sitting under 
gravity stress. The slope of the compliance curve must, 
therefore, be made smaller as per the following example: 

0LD CARD . ... 

6 0 • -5* 800 • 40* 1000* 80* 1200* 140* 2000* 

NEW CARD 

6 0 • -5. 9 0 0 • " 40* ~ i 000 • ~50'* HOD* i 4 O'* 14 00* 


♦Additional information is available in Appendix VI. 


190* 


190 "* 
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Note that since a nominal pressure in the prone legs of 
80 mm of mercury is to be maintained, the curve must pass 
through the initial leg volume of 80 mm of mercury. 

B. Change in Volume of Kidney 

Compliance changes must always accompany volume changes 
unless fluid loss or gain is to be simulated. Suppose we wish 
to keep the model basically the same, but reduce kidney volume 
by 50%. The following changes will accomplish the desired 
result : 

9LD CARDS 

4 0* "5* 45q* 40* 500* 8o« 650* 140* 1000* 130* 

24 500 

NEW CARDS 

4 0* *5. 225* 40* 250* So* 3?5* 140* 500* 190* 

24 250* 


C. Fluid Loss in System 

The simplest way to simulate fluid loss is to reduce 
initial volume without changing compliance. While the blood 
can be removed from any organ, or by taking part of the blood 
out of all organs, a simple way is to remove the total amount 
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to be depleted from the vena cava. The following illustrates 
a 150m£ blood loss: 


6LD CARO 
27 350 • 

NEW CARD 
27 200 • 


It should be remarked here that since the model does 
not consider intravascular or extravascular fluid shifts, 
removal of 100 to 15 0 mil of blood from the model is about 
equivalent to removal of h liter of blood from a human. 

D. Increase Flow to Legs to Simulate Exercise 

Since F = ?2 , we can double flow by halving 

R 

resistance. The following change will double the flow 
through the legs. There will, however, be a corresponding 
rise in cardiovascular output and an increase in heart rate 
automatically to accomplish the desired change. 

9LD CARD 
76 6*00.0014 

NEW CARD 
76 3.00.0014 
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E. Change in Cardiovascular Output without Heart Rate Change 
or vice versa 

These changes require a considerable fine tuning of 
all organs from the vena cava through the last section of 
the aorta. Resistances and compliances must be changed 
together throughout the system, accompanied by inertial 
changes to maintain realistic flows and especially to 
maintain the correct amount of back flow in the aorta. No 
attempt will be given here to explain these changes. They 
are best accomplished by the trial and error method and had 
best be made in small amounts. Often a 10% change in value 
of compliance can be easily seen on a plot, and often a 50% 
change will cause a mis-balanced system that will no 
longer operate. 

Consulting any good medical physiology textbook will 
yield considerable insight into the behavior of this model, 
both in normal and "pathological" modes of operation. 
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APPENDIX I 


DERIVATION OF AORTIC PRESSURE-VOLUME EQUATION 
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APPENDIX I 


DERIVATION OF AORTIC PRES SURE- VOLUME EQUATION 


1) T 


E — + AS. dA r 

r r dt 
o o 


Pr 

6 


2 ) 


Ar = r-r 

o 


(— ) r 
r o 


E + ( 


Re dr 

r ' dt 
o 


Pr 

6 


3) (f) - (K«)± + £ 

o o 


= P 


A = frr 


dr _ _1_ dA 
dt 2Trr dt 


4) 



■nh 

< E6 )Ai 


6 Re 1 dA 

+ ( r ) 2iTr 2 dt 


P 


5) 


(f)a-(f)S + 

o 


6 Re 1 dA 

( 2r ''A dt 
o 


P 


By Taylor Series Expansion 


6 ) 


i-(¥>* -"4 (4- 


i) 
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DERIVATION OF AORTIC PRESSURE-VOLUME EQUATION (CONTINUED) 


7) 


5 E 6 

( 4 r A 
o o 


)(A-A 0 ) + C 


6 Re 1 dA 

2r ' A dt 
o 


P 


or 


E 

s 


5 E 5 
4 r A 


o o 


E 


d 


(SRe 

2r 

o 


1 dA 

8) E (A-A ) + E T dt = P 
sod 


But, V = A*X where X is length of each section of the aorta. 
Thus , 


9) 


E 

S X 


1 dV 
+ E d V dt 


P 


By adding the same gravity term as in basic pressure, equation 4) 
yields : 


10 ) 


P = E 


lYzloi 

X 


+ G + 


Id 

V 


dV 

dt 
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APPENDIX II 

DESCRIPTION OF VARIABLES USED IN COMPUTER PROGRAM 
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APPENDIX II 

DESCRIPTION OF VARIABLES USED IN COMPUTER PROGRAM 


VARIABLE"/ UNITS/ 

NORMAL range/ 

DESCRIPTI0N 

^ADDF 

CM##3/SEC 

•100 TO 2000 

FLOW TEMP STORAGE „ 

HEART FLOW (CARDIAC OUTPUT) TEMP STORAGE 

REAL REPRESENTATION OF I I __ 

bbdy angle with convention o • horizontal 
SURFACE AREA OF rODY NODES 

AFL0 

* All 

CM**3/SEC 

INTEGERS 

5q TO 400 
1 TO 1200 

angle 

aold 

degrees 

CM*#2 

0 To 90 
10 TO 1000 

APRS 
• ARAD 

MM HG 
RADIANS 

90 TO 110 
0 TO 6*28 

mean blood pressure temp storage 

BODY angle WITH CONVENTION 0 » HORIZONTAL 

Compliance table ( see text for detail's ) 

CURRENT TIME 

DERIVATIVE OF AREA WITH RESPECT TO TIME 
TIME F0R ONE COMPUTE CYCLE 

’ CqMP 

' CTIME 

VARIOUS 

SECONDS 

•100 T0 2000 
0 TO 36000 

dadt 

delt 

CM**2/SEC 

SECONDS 

•100 T0 100 ' 
O.OOl TO 0*01 

DHRATE 
i DIAS 

beats/min 

MM HG 

•20 TO 20 
5o TO 100 

DELTA HEART RATE FOR CONTROL OF HEART 
DIASTOLIC BLOOD PRESSURE TEMP STORAGE 

DIST 
. - D0T6 

CM 

SECONDS 

•200 TO 200 
•000b TO *005 

body node distance from heart 

ONE HALF TIME FOR ONE COMPUTE CYCLE 

DTIME 

DVDT 

MILLISECONDS 

CM»*3/SEC 

0 TO 1200 
•1000 TO 1000 

TIME WITHIN BEAT 

DERIVATIVE OF VOLUME WITH RESPECT TO TIME 

EDP 

. ERTIA 

MMHG*SEC*"*2/CM**3 *4 T0 4 
MMHG*SEC»*2/CM**3 •0005-* 

“INERTIAL VALUE FOR AORTIC WALLS" 

002 INERTIAL VALUE FOR BLOOD 

ES p 
F 

“NOT USED IN THIS FORM 
CM**3/SEC -1OQT02OOO 

static inertial value of aortic walls 

OUTPUT BUFFER OF VALUES OF FOLD 
COMPUTE STEPS PER HEART BEAT 
FLOW BETWEEN BODY NODES 

» FI 

FOLD 

INTEGERS 

CM#*3/gEC 

10 TO 240 
•100 T Q 2000 

FpRV 

■ ^^FTIME 

CM**3/SEC 

seconds 

•100 TO 2000 
0 TO 36000 

PREVIOUS FLOW BETWEEN BODy NODES 
CURRENT time temp storage 

^■grav 
» W G(1) 

G 

MM HG 

0 TO 5 
•50 T0 50 

GRAVITY FIELD 

LOWER BODY PRESSURE. OUTPUT VECTOR 

' G(2) 

' G ( 3 ) 

MM HG 
MM HG 

100 TO 200 
50 T 0 100 

SYSTOLIC BLOOD PRESSURE OUTPUT VECTOR 
DIArTbLIC BLBOD PRESSURE OUTPUT VECTOR 

G ( 4 ) 
G ( 5 > 

MM HG 
beats/min 

90 TO HO 
0 TO 180 

MEAN b^OOD PRESSURE OUTPUT VECTOR 

heart rate output vector 

G ( 6 ) 
G ( 7 ) 

DEGREES 

CM##3/SEC 

0 TO 90 
50 TO 400 

BODY ANGLE OUTPUT VECTOR 

HEART FLOW (Cardiac OUTPUT) OUTPUT VECTOR 

G ( 8) 
HAFLO 

G 

CM**3/SEC 

0 TO 5 
50 TO 400 

GRAVITY FIELD OUTPUT VECTOR 
HEART FLOW (CARDIAC OUTPUT) 

HAPRS 

HDIA S 

MM HG 
MM HG 

90 TO 110 
50 TO 100 

MEAN BLOOD PRESSURE 
DIASTOLIC BLOOD PRESSURE 

HI qH 
HMAX 

CM 

BEATS/MIN 

•200 TO 200 
0 TO 180 

BODY NODE HEIqHT FROM HEART 
MAXIMUM HEART RATE 

HMIN 

HOLD 

BEATS/MIN 

RE A l 

0 TO 180 
ANY 

MINIMUM HEART RATE 

INPUT VECTOR FOR ALL INPUT CONSTANTS 

HRATE 

HSYST 

beats/min 
MM HG 

0 TO 180 
100 TO 200 

HEART HATE 

SYSTOLIC BLOOD PRESSURE 

HTIME 

I 

SECONDS 

INTEGER 

0.3 TO 1.2 
ANY 

t t he ffJR 0 ne heart beat 
LOOP COUNTER 

IC0MP 
I EC 

1/SEC 

MILLISECONDS 

100 TO 1600 
0 TO 1200 

COMPUTE FREQUENCY INTEGER 
TIME WITHIN BEAT 

IETMS 

IFIRST 

IH0LD 

MILLISECONDS 

INTEGERS 

INTEGER 

INTEGERS 

0 TO 1200 

1 T0 10000 
"ANY 

1 TO 1200 

beat time 

BUFFER COUNTER 

NUMBER OF PHYSICAL UNIT OF STORAGE TAPE 
COUNTER FOR COMPUTE LOOP 

^Fin 

INDEX 

INTEGER 

UNITLESS 

ANY 

N0NE 

NUMBER OF PHYSICAL UNIT OF CARD READER 
NOT currently used 

I NDX 

integer 

ANY 

INTERMEDIATE VARIABLE 




... . - - - — 
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DESCRIPTION OF VARIABLES USED IN COMPUTER PROGRAM (CONTINUED) 


_ VARIABLE/ UNITS/ 

NORMAL RANGE/ 

^Biotest 

integer 

any 

iout 

INTEGER 

ANY 

ioutf 

1/SEC 

50 TO 200 

IPRINT 

INTEGERS 

0 OR l 

IRUN 

INTEGERS 

_any 

I SKIP 

INTEGERS 

1 TO 16 

ISTEP 

INTEGERS 

10 TO 240 

ITEST 

"integer 

0 TO 99 

IT IME 

MILLISECONDS 

3q0 TO 1200 

J 

integer 

ANY 

K 

integer 

ANY 

LIEC 

MILLISECONDS 

0 TO 1200 

LOOP 

integers 

200 TO 800 

N 

INTEGER 

"ANY 

NT IME 

MILLISECONDS 

300 TO 1200 

P 

MM HQ 

-5 T'0 200 

PADJ 

MM HG 

•5 TO 200 

• PLBNP 

MM HG 

*50 TO 50 

) POLO 

MM HG 

*5 TO 200 

PSUBR 

MM "HQ 

■5 

‘ RE 

MMHQ#SEC/CM**3 *0001 TO 10 

| SANG 

UN IT LESS 

-1.0 T 9 1.0 

• SINC 

1/SEC 

100 TO 1600 

SYST 

MM HG 

100 TO 200 

> ^^TACC 

SECONDS 

1 TO 36000 

’ ^Atemp 

1/SEC 

50 TO 200 

1 ^^TIME 

seconds 

1 TO 6000 

> TlNC 

SEC 

0.1 TO 10 

■ TNOW 

seconds 

1 TO 36000 

1 TOUT 

seconds 

1 TO 36000 

< TVOL 

CM**3 

3.000_T0_6p00 

i UPDATE 

CM**3 

*2* TO 2* 

t V 

CM**3 

10 T0 2000 

i VC 

CM**3 

10 TO 2000 

1 VOLD 

CM**3 

10 TO 2000 

- VPRV 

CM**3 

10 TO 2000 


D ESC RIPTION 

INTERMEDIATE VARIABLE _ 

NUMBER OF PHYSICAL UNIT OF PRINTER 

OUTPUT FREQUENCY I NTEG ER 

print option switch# o for no print 

NUMBER OF CURRE NT RUN _ 

ICOMp/IOUTF 

COMPUTE STEPS PER HEART BEAT 

INTERMEDIATE STORAGE OF INpgT CARD NUMBER 

INTEGER TIME FOR ONE HEART BEAT TEMP 

LOOP COUNfER 

L0OP_ COUNTER 

TIME OF END OF BEAT 

COMPUTE STEPS PER OUTPUT BUFFER 

TEMP "STORAGE FOR INPUT 

INTEGER TIME F8R ONE HEART BEAT 

PRESSURE AT NODES# "OUTPUT BUFFER 

PRESSURE ADJUSTMENT FOR ORGAN HEIGHTS 

LOWER BODY PRESSURE# ^"VALUES HAVE -SIGN 

PRESSURE AT NODES 

PRESSURE OF COLLAPSED WALL 

RESISTANCE TO FLOW BETWEEN NODES 

SINE OF BODY ANGLE 

COMPUTE FREQUENCY 

SYSTOLIC BLOOD PRESSURE TEMP STORAGE 

LENGTH OF CURRENT RUN 

OUTPUT FREQUENCY 

LENGTH OF current TIME segment 

time for filling one output buffer 

CURRENT TIME WITHIN SEGMENT 

OUTPUT VALUE" OF "CURRENT TIME 

TOTAL VOLUME OF BLOOD 

ERROR IN BLOOD VOLUME BEFORE UPDATE 

OUTPUT 3UFFER 0F VALUES OF vOLD 

variable USED for VOLUME CHECK 

VOLUME OF BODY NODES 

VOLUME OF BODY NODES AT LAST STEP 
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APPENDIX III 
FLOWCHARTS 
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GENERAL PROGRAM FLOWCHART 
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INNER COMPUTATION LOOP FLOWCHART 
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APPENDIX IV 


COMPLIANCE CARDS 
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APPENDIX IV 


COMPLIANCE CARDS 


The compliance cards 1-16 are read with a format of I2,10F7.2. 
The 12 field is the same two digit indicator field as other control 
cards. Numbers 1-15 are compliance for the 15 organs. Number 16 is 
a heart rate control table. 


The other 10 fields should be considered as five pairs of values. 
The first of each pair {odd values) is in all cases the independent 
variable. The second of each pair (even values) is the dependent 
variable. The table is always entered with an independent variable. 

The dependent variable is computed by linear interpolation and returned 
to the program for computation. 


The units for these variables are in three classes : 


Indicator Field 


Independent Variable 


Dependent Variable 


16 


Pressure in mm Hg at 
baroreceptor 


Beats per minute to be 
added to heart rate 


8,11,12 


Time in milliseconds 
within a heart beat 


Compliance C in 


CM 

mm Hg 


All other 
fields 


Volume V of organ in Pressure P in mm Hg 

cm 3 


Note that for the heart chambers 8, 11, and 12, since P is a function 
of both volume and time, a more complicated computation is necessary. 

In the other chambers, since P is a function of volume only, we save 
considerable time at no loss of generality by returning P instead of C. 
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COMP LI AN CE CARDS (CONT IN UE D ) 


An example of use follows. The following is the image of the 
compliance card for the legs : 


6 0* -5 * 800* 40* 1000* 80 1200* 140* 2000* 190t 


For a value of 1000 cm 3 , the pressure returned will be 80 mm Hg . 
For 1200 cm 3 , the returned pressure will be 140 mm Hg. If the input 
volume were 1020 cm 3 , then by linear interpolation the pressure returns 
will be 86 mm Hg. 

The current compliance table is given below for reference. 


1 

0* 

■5 * 

45q • 

40* 

500* 

80* 

650* 

140. 

1000* 

190 

2 

o* 

"5 • 

450* 

40 • 

500 • 

so* 

650* 

140' 

1000* 

190 

3 

0* 

-5. 

450* 

40 • 

500* 

80* 

650* 

140 

1000 *. 

190 

4 

0* 

-5* 

45q • 

40* 

500* 

80* 

650 * 

140* 

1000 

190 

5 

0* 

*»5 « 

45o» 

40* 

500* 

80* 

650* 

140* 

1000* 

190 

6 

O 

•5 • 

800* 

40 • 

1000* 

80* 

1200* 

14Q. 

2000* 

190 

7 

0* 

*5 • 

200* 

5 * 

475. 

10* 

750* 

15. 

2000*. 

90 

8 

O 

20 • 

50 * 

2.7 

170* 

2.2 

220* 

100* 

500*. 

■ 90 

9 

0* 

-5. 

45Q. 

0* 

500* 

40* 

650 » 

100* 

1000. 

150 

10 

0* 

*5 * 

450 

0* 

500* 

4 0. 

650. 

100* 

1000* 

150 

11 

0* 

4.5 

6Q * 

3.5 

90* 

500* 

130* 

13. 

1200*. 

8 

12 

100* 

24* 

1 5 0 » 

1.4 

250* 

1.1 

330. 

48. 

700* 

?4 

13 

0* 

*5 * 

20* 

65* 

30* 

100* 

50* 

ISO* 

200* 

200 

14 

0* 

*5 • 

20 • 

65* 

30* 

100* 

50. 

iso. 

200* 

200 

15 

0* 

■ 5 * 

20* 

65* 

30* 

100* 

50* 

iso. 

200* 

200 

16 

85 • 

25. 

97* 

1* 

99* 

•01 

101* 

“ * 01 

121* 

-?0 
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APPENDIX V 


DESCRIPTION OF VARIABLES USED IN TEXT 
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APPENDIX V 


DESCRIPTION OF VARIABLES USED IN TEXT 
SYMBOL DESCRIPTION 


P 

P 

P 

V 

V 


pressure 

pressure for collapsed 
pressure at organ n 
volume 

volume at time o 


o 


F 

n 

F 

C 

t 

R' 

R 

X 

A 


r 

E s 

E d 

L' 

L 

P 

n 

6 

T 

E 

R e 

G 

t 


flow rate 

flow rate at step n 
compliance 
t ime 

resistance/unit length 

resistance 

length 

area 

area at time o 
radius 

static elasticity 

dynamic elasticity 

density • length 

inertial constant for 

density 

viscosity 

wall thickness 

wall tension 

Young's Modulus 

wall drag coefficient 

gravity 

time 


wall 


blood 



i-i \o> 


DESCRIPTION OF VARIABLES USED IN TEXT t CONTINUED) 


Assumed Relationships 


<$n 
Trr 4 

o 

ir r 2 

Ar Al. 

T = E— + R d r 

r dt 


R' 

L' 
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USER'S GUIDE 

for 

LBNP PROGRAM AND THE ASSOCIATED DISPLAY PROGRAM 
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USER'S GUIDE 
for 

LBNP PROGRAM AND THE ASSOCIATED DISPLAY PROGRAM 

GENERAL 

This program is a model of the human cardiovascular system. 
The model consists of 15 body compartments or organs and the 21 
associated flow paths between these organs (Figure 1) . The body 
parameters in the sample deck supplied with the program are for 
a "typical" male human. That is to say it functions as a male 
of average height and weight. These parameters have been tested 
under various conditions such as positive and negative lower 
body pressure, various G loading between -2 and +10 G's, and 
various postures such as lying, standing, and sitting. In all 
cases, parameters such as heart rate, flow, cardiovascular 
output, pressures, and volumes are about what one would expect 
with a conditioned man. 

The basic formulae for the model are: 



where V is volume, P is pressure, F is flow, R is resistance, 
C is compliance, G is gravity, and L is inertia. 
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Figure 1 

MODEL CONFIGURATION 
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Flow to any organ may be altered by changes in resistance 
between the organs with a smaller resistance yielding a larger 
flow. In this way exercise may be simulated. If we reduce the 
resistance to the arms to one-half we will get approximately 
double the normal flow to the arms along with the associated 
increase of both heart rate and cardiovascular output. 

Equation 2) is expedited by use of a volume pressure function 
where for a given volume, a value of P is returned. The function 
for the legs of a man with good conditioning and 1000m& blood 
volume in the legs may look like Table la. 

Conditioned 

Volume Pressure 

0 -5 

800 40 

1000 80 

1200 140 

2000 190 

Table la Table lb 

Thus, for the unconditioned man a larger volume change is 
necessary to obtain a given pressure. Thus, if the unconditioned 
man is subjected to standing under any G forces, more blood will 
tend to flow to his legs before the model stabilizes. In most 
cases the model will stabilize in 20 seconds after minor changes, 

ZB- 
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but may take 60 or more seconds if large changes are made. 

The input of physiologically unreasonable data constants 
will, in general, create conditions that will cause the program 
to fail. In general, these same conditions if applied to the 
human body will also cause it to fail. Examples of this are 
given, but this list is not meant to be exhaustive, namely: 

a) Changes in initial volume of an organ must be 
accompanied by changes of the compliance table. 

b) Changes of resistance and inertia affect flow 
rates and volume changes. Thus, changes should 
be relatively small; and in many cases, they 
should be accompanied by changes in compliances . 

c) Rapid discrete . changes in body angle will result 
in immediate changes of pressure, and thus sharp 
changes in flow will occur. 

d) There are no default options for volumes, pressures, 
compliances and resistances. Failure to input 
these constants leads to failure of the model. 

e) Because the computation is in discrete space instead 
of a continuum, the compute frequency is critical. 
Generally, at least 200 steps per second should be 
considered the minimum compute frequency. One 
should bear in mind that running time is proportional 
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to compute frequency. 


The general input and output for the program is very 
flexible. For example: 

a) Parameter cards may be read in any order. 

b) Parameter cards can be duplicated. Only the 
last card read will be used. 

c) A parameter, once read, will remain unchanged 
until another parameter of the same type is 
read. 

d) Any of the computed parameters may be output 
graphically, with from 1 to 20 of these 
parameters displayed in a single run. (Up 
to 40 parameters are allowed if a 22 inch 
plotter is used) 

e) The display program is separate from the compute 
program s'o that any number of display runs may 
be made from a single compute run. 

f) Both time and value scaling of the graphical 
output is flexible and independent of the 
output frequency. 

g) Output parameters may be displayed in any order. 



DECK SETUP FOR LBNP PROGRAM 


Input cards are as follows: (Format (12, 10F7.2)) 


TABLE II 


COMPARTMENTS 

1 BRAIN 

7 

VENA CAVA 


13 

AOR- 

2 

ARMS 

8 

RIGHT HEART 


1A 

A0R 

3 

Trunk 

9 

RIGHT LUNG 


15 

A0R 

4 

kidney 

10 

LERT LUNG 




5 

GUT 

11 

LEFT ATRIUM 




6 

LEGS 

12 

LEFT VENTRICLE 




FLOW RELATIONSHIP 

BETWEEN 

compartments, v 

DENOTES 

i VALVE* 

1 

1 "> 7 

V 

8 8 -> 9 

V 15 

8 

•»1D 

2 

2 •> 7 

V 

9 9 *>11 

16 

1 5 

•> 1 

3 

3 -> 7 

V 

10 10 ->11 

17 

15 

- > 2 

4 

4 •> 7 

V 

11 11 ->12 

18 

15 

-> 3 

5 

5 *> 7 

V 

12 12 ->13 

19 

15 

■•> 4 

V 6 

6 -> 7 


13 13 ->1A 

20 

15 

-> 5 

V 7 

7 •> 8 


14 14 ->1 5 

21 

15 

-> 6 


INPUT CARDS 

1-15 

16 

21 - 35 
41 - 5b 
60 

61 - 66 

70 

71 - 91 

92 

93 

94 
9b 

96 

97 

98 

99 


COMPLIANCE CARDS 

MEAN PRESSURE*DELTA HEARTRATE TABLE* 

INITIAL VOLUME CARDS 
INITIAL PRESSURE CARDS 
GRAVITY IN G'S* 

DISTANCE FROM HEART T9 9RGANS 1-6 IN CM. 
INERTIA C9NSTANT F0R A0RTA. 

RESISTANCE AND INERTIA BETWEEN N0DES 
INITIAL HEART RATE IN BEATS PER MINUTE. 
L9WER BBDY NEGATIVE PRESSURE IN MM HG* 

ANGLE IN DEGREES. 0 LYING# 90 STANDING* 
TIME IN SECONDS FSR RUN WITH ABOVE PAR. 
COMPUTE ANq OUTPUT FREQUENCY. 

PRINT OPTION. 0 NOpRINT# 1 PRINT FIRST BUF 
END 9F RUN. 

END 9F JOB. 


The program will read cards until a 95 card (Time) , a 98 
card (End of Run) or 99 card (End of Job) is read. A convenient 
way of using the program is as follows: 
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Figure 2 

DECK SETUP FOR LBNP PROGRAM 


99 End of Job Card 


98 End of Run Card 


95 Time Card 


DATA CARDS 


95 Time Card 


DATA CARDS that are user supplied; 
i.e., modifications of the 
existing model. 


DATA CARDS supplied with program 

up to but not including 

the first 95 Card. 


LBNP Program 


Multiple or stacked runs may be run by inserting more 
Data cards and Time cards beyond the first 98 End of Run 
card. Each run must be ended with a 98 End of Run card. 

If used in this way, no parameter cards will be missing 
If a change to the model is considered permanent then it can 
be put in the place of its counterpart in the Data cards 
supplied with the program. 

The program will compute and produce output for the 
number of seconds specified on the Time card. This time 
should be whole seconds if an output frequency of 200 or 
more is used, even seconds if an output frequency of 100 
or more is used, 4 times an integer number of seconds if 50 
is used, etc. After the specified number of seconds, the 
program will read new parameter cards until a new Time card 
is read. Only changes need be added, since all unchanged 
parameters will remain fixed. The LBNP will create a data 
tape to be read by the DISPLAY Program. Compute frequency 
should be a non-negative integer power of 2 multiple of the 
specified output frequency; i.e., if output frequency is 
100, then compute frequency should be 100, 200, 400, 800, 
etc. 
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DECK SETUP FOR DISPLAY PROGRAM (Figure 3) 

Item 

1. DISPLAY Program 

2. Scale Cards supplied with Display Program in order 
received. 

3a. One card specifying number of output graphs, maximum 
graph height, and seconds per display inch in 
FORMAT (14 , 2F4 . 0) . 

NOTE: 

Number of graphs must be between 1 and 20. 

(up to 40 are allowed for 22 inch plotter) 

Height must be hr 1 or 2 inches. 

(up to 4" allowable with 22 inch plotter) 

Seconds per inch can be any real number. 

3b. One card specifying which parameters are to be displayed, 
FORMAT (2014); the parameter codes are listed in Table II. 

For 22 inch plotters, if more than 20 parameters are to be 
displayed, then two of the above cards must follow each 
card specifying the number of graphs . 
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TABLE III 


DISPLAY C6DES F8R PARAMETERS 


101-121 

FL8wS 1-21 

201-215 

VOLUMES 1-15 

301-315 

PRESSURES 1-15 

401 

IBNP (L9WER BODY PRESSURE) 

402 

SYSTQLIC PRESSURE 

403 

01AST9LIC PRESSURE 

404 

MEAN BL99D PRESSURE 

405 

heart rate 

406 

B90Y angle 

407 

CARO I 9 VASCULAR OUTPUT 

408 

GRAVITY 

411 

SYST9LIC* MEAN AND DIASTOLIC PRESSURE 



LO 




RUNNING TIME 


On the XDS Sigma 5 computer with a compute frequency of 
200 steps/second and an output frequency of 50 steps/second, 
the execution time of the LBNP Program is 3.373 minutes for 
the enclosed experiment which lasted 56 seconds. This amounts 
to about 3*s times as much computer time as the desired model 
time. The execution time of the Display Program to generate 
the enclosed plot is 1.280 minutes, or nearly real-time. The 
digital plotter required 4.565 minutes to plot this display. 
The same display at % inch per second would require about 
twice this amount of compute and plot time. (See Appendix) 

Remark on Stepsize : As seen in the previous section, the 

Basic Flow equation (5) becomes a first order, first degree 
differential equation (7) when inertia is added. In the 
descrete approximation to (7) the corresponding difference 
equation 

P - P - G - L 
p = m n A t 


R 

m,n 


is sensitive to stepsize changes. Larger than real life 
inertial values and greater than real life flow rate changes 
create an unstable situation for large stepsizes. In general, 
a high-frequency ringing of flow, especially in the aorta, 
can be stopped only by using a smaller stepsize (higher compute 
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frequency) or by reducing the blood inertial values. A 
compute frequency of 200 is usually sufficient for a stable 
condition with any realistic inertial values or flow rates. 

The DISPLAY Program reads the tape created by the LBNP 
Program and produces the specified graphs . There is no 
provision for multi _ runs from the same data tape; however, 
the program may be run any number of times using different 
plot parameters with the same data tape. 

For information concerning this model that is not in this 
user's guide, see "A Cardiovascular System Model for Lower- 
Body Negative Pressure Response" by Mitchell, B.A. and Giese, 
R.P., included with the Quarterly Report No. 5, NASA Contract 
No. NAS 9-11119. 
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APPENDIX 


UNI VAC 1108 DECK SETUP 
for 

RUNNING BOTH PROGRAMS TOGETHER 
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UNIVAC 1108 Deck Setup for Running both Programs Together 


Z RUN B02380#D36i243,8636N,Q5o9,C* 10/5 MAULDIN 

N MSG FILE READ TAPE 3 FH432 0 FSTrAN 1 

S ASG A*58500 
SR ASG C-63402 
ASG F-CALPLT 
XQT CUR 
TRW A# Ci F 
IN C 
TRI C 
F0R A 1 # A 1 

(LBNP PR0GRAM) 

F0R AH/AH 


(DISPLAY PROGRAM) 


XQT A 1 


1 

0 • 

2 

0* 

3 

0* 

4 

0* 

5 

0* 

6 

0* 

7 

0* 

8 

0* 

i 9 

0* 

f 10 

0* 

11 

0* 

12 

100* 

13 

0* 

14 

0* 

15 

0* 

16 

■85 • 

21 

500* 

22 

500* 

23 

500* 

24 

500* 

25 

500 • 

26 

1000* 

27 

350* 

28 

250* 

29 

500* 

30 

500 • 

31 

120* 

32 

120* 

33 

30 • 

34 

30* 

35 

30* 

41 

80* 

42 

80* 

k 43 

80 • 

1 44 

80* 

f 45 

80* 

46 

80 • 

47 

5* 

48 

100* 

49 

50* 


■5 • 

450 • 

40 • 

500* 

Bo* 

650 • 

140* 

1000* 

190* 

•5 • 

450* 

40* 

500* 

8o • 

650 • 

140* 

1000* 

1?0« 

•5 • 

45o» 

4Q • 

500* 

So* 

650 • 

140* 

1000* 

1?0« 

• 5 * 

45o.» 

40* 

500 • 

8o . 

650 • 

140* 

1000* 

190* 

-5. 

450* 

40* 

500* 

80* 

650 * 

140» 

1000* 

190* 

•5 * 

800* 

40* 

1000* 

«o* 

1200* 

140* 

2000. 

190* 

•5 • 

200 • 

5. 

475* 

10* 

750* 

15. 

2000* 

90* 

20 • 

50» 

2.7 

170* 

2.2 

220* 

100* 

500* 

20* 

-5. 

45o» 

0* 

500* 

40. 

650* 

100. 

1000* 

150* 

•5 • 

450 * 

0» 

500 • 

40* 

650 * 

100. 

1000. 

150* 

4.5 

60* 

3.5 

90 • 

500. 

130* 

13. 

1200. 

8. 

24. 

150. 

1.4 

250* 

1.1 

330* 

48. 

700. 

?4 . 

"5 . 

20* 

65. 

30* 

100* 

50* 

iso* 

200* 

290* 

• 5 . 

20* 

65 * 

30. 

100. 

50* 

iso. 

200* 

290* 

•5 . 

20* 

65* 

30* 

100* 

50* 

iso. 

200* 

290* 

25. 

97- 

1 • 

99. 

•01 

1 0 1 • 

•*01 

121* 

-20* 



50 

50* 

51 

100* 

52 

100* 

53 

100* 

54 

100* 

55 

100. 

61 

35. 

62 

20. 

63 

-15. 

64 

*22*5 

65 

*30* 

66 

• 70 * 

70 2.5 

71 

4.00*0014 

72 

6.00*0014 

73 

6.00*0014 

74 

6*00.0014 

75 

6*00*0014 

76 

6.00*0014 

77 

0*02*0014 

78 

0*30.0014 

79 

0*45*0014 

80 

0.45.0014 

81 

0*003*0007 

82 

0*005*0014 

83 

0*007.0014 

84 

0*010*0014 

8b 

0*30.0014 

86 

1.00*0014 

87 

1*50*0014 

88 

1.50*0014 

89 

1*50*0014 

90 

1.50*0014 

91 

1*50*0014 

92 

80* 

93 

o. 

94 

o. 

96 

200. 50. 

95 16. 

93 

-30* 

95 

8. 

93 

-40* 

95 

8* 

93 

•50* 

95 

8* 

93 

0* 

95 16. 

98 


99 


XQT AH 

101 

0 *02 

102 

0 *05 

103 

0 .05 

104 

0 *05 

105 

0 *05 

106 

0 .02 

107 

0*002 

108 

0*002 

109 

0 *01 

110 

0 *01 



111 

0* 

002 

113 

0* 

001 

113 

0* 

001 

114 

ot 

001 

115 

0» 

001 

116 

0 

• 01 

117 

0 

• 01 

118 

0 

• 01 

119 

0 

• 01 

120 

0 

• 01 

121- 

100* 

005 

201 

100* 

002 

202 

100* 

002 

203 

100* 

002 

204 

100* 

002 

205 

100* 

002 

206 

0* 

001 

207 

0» 

001 

208 

0* 

002 

209 

400* 

005 

210 

400* 

005 

211 

0* 

005 

212 

0* 

002 

213 

0 

• 01 

214 

c 

•01 

215 

0 

• 01 

301 

0 

• 01 

302 

0 

• 01 

303 

0 

• 01 

£04 

0 

• 01 

?305 

0 

• 01 

306 

0« 

005 

307 

0 

.05 

308 

0* 

005 

309 

0 

• 01 

310 

0 

• 01 

311 

0 

• 02 

312 

0* 

005 

313 

0* 

005 

314 

0* 

005 

315 

0» 

005 

401 

0- 

■•01 

402 

0* 

005 

403 

0* 

005 

404 

90 

• 05 

40b 

50 

•01 

406 

0 

•01 

4Q7 

0 

•01 

408 

0 

• 10 

409 

0 

• 01 

410 

0 

• 01 

411 

0* 

005 

20 

0*5 

10 

4q5 

4q7 

401 


101 


116 ?oi 301 106 1 PI 206 3q6 20? 30? 5 0* ?09 P13 313 314 315 41 


I TRI 
E PMO 
EOF 


A# J 
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PRES 


P 15. 
P 14. 
P 13. 

V 13. 

V 9. 

V 8. 

P 7. 

V 7. 

P 6. 

V 6. 

F 21. 
F 6. 

P 1. 

V 1. 

F 16. 
F 1. 
LBNP 
CARD 
HR 


200 . 


0. 

200 ‘ 

0 . 

2 °°' 

0. 

200 . 

0. 

100 . 

0 . 

600 . 

400 . 

500 . 

vvvVYY'vrv'VVVVVHA/VVYVVVvvvvvvvvvvvvwvwwvwwvvw/vvvvvvvvvv^ 

0. 

20 . 


0. 

1000. 



1000. 


0. 

100 ‘ 

- 100 . 

50 . 


0 . 

100 . 


0 . 

600 . 


100 . 



SO. 


0. 

- 100 . 


0. 

100 . 


0. 

150 . 
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